Does the NJL chiral phase transition affect the elliptic flow of a fluid at fixed rj/s? 
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We have derived and solved numerically the Boltzmann-Vlasov transport equations that includes 
both two-body collisions and the chiral phase transition by mean of NJL-field dynamics. The scope 
is to understand if the field dynamics supply new genuine effects on the build-up of the elliptic 
flow i>2 , a measure of the asymmetry in the momentum space, and in particular if it can affect the 
relation between V2 and the shear viscosity to entropy ratio n/s. Solving the transport equation with 
a constant cross section for the condition of Au + Au collisions at y^ijvjv = 200 AGeV it is shown 
a sizable suppression of V2 due to the attractive nature of the field dynamics that generates the 
constituent mass. However the key result is that if rj/s of the system is kept fixed by an appropriate 
local renormalization of the cross section the V2 does not depend on the details of the collisional 
and/or field dynamics and in particular it is not affected significantly by the chiral phase transition. 

PACS numbers: 25.75.-q , 25.75.Ld, 12.38.Mh, 24.85.+p 



Introduction - The ultra-relativistic heavy-ion colli- 
sions at high energy ^Js^n ~ 200 AGeV represent the 
main tool to study the formation and the properties of 
the quark-gluon plasma at high temperature. The RHIC 
program at BNL has shown that the azimuthal asymme- 
try in momentum space, namely the elliptic flow 1)2, is 
the largest ever seen in HIC suggesting that an almost 
perfect fluid with a very small shear viscosity to entropy 
density ratio, rj/s, has been created From simple 

quantum mechanical considerations [4| as well as from 
the study of supersymmetric Yang-Mills theory in the 
infinite coupling limit a lower bound for rj/s of about 
~ 10 _1 is predicted [5]. Such a value is much lower than 
any other known fluid and in particular smallerjman the 
one of water and even than the superfluid He [f|. 

First developments of relativistic viscous hydrodynam- 
ics @, [H as well as parton cascade models [9MTl1] indicate 
that even a small rj/s ~ 0.1 — 0.2 affects significantly the 
strength of V2(pt) especially at pr > 1 GeV. Therefore 
it has become mandatory to determine the value of rj/s 
of the plasma created at RHIC through the study of the 
relation between rj/s and V2 [12| . However viscous cor- 
rections to ideal hydrodynamics are indeed large and a 
simple relativistic extension of first order Navier Stokes 
equations is affected by causality and stability patholo- 
gies [H, [l4j]. It is therefore necessary to go to second 
order gradient expansion, and in particular the Israel- 
Stewart theory has been implemented to simulate the 
RHIC collisions providing an upper bound for rj/s < 0.4 
[lf|. Such an approach, apart from the limitation to 
2+1D simulations, has the more fundamental problem 
that it is based on a gradient expansion at second or- 
der that is not complete [l3j]. Furthermore it cannot be 
sufficient to describe correctly the dynamics of a fluid 
with large 77/s as the one in the cross-over region and/or 
hadronic phase which at least at RHIC still gives a non 
negligible contribution to V2 [H| that affects the deter- 



mination of the rj/s itself YV\. 

A relativistic transport approach has the advantage 
to be a 3+ ID approach not based on a gradient ex- 
pansion that is valid also for large viscosity and for out 
of equilibrium momentum distribution allowing a reli- 
able description also of the intermediate pt range where 
the important properties of quark number scaling (QNS) 
of v 2 (pt) have been observed (l8| . In this px region 
viscous hydrodynamics breaks its validity because the 
relative deviation of the equilibrium distribution func- 
tion 5f I f eq increases with p\ becoming large already at 
p T > 3T~ lGeV [H. 

In this perspective transport approaches at cascade 
level have already been developed |9l-fIiT |20| . but they 
miss any effect of the field interactions responsible for 
the chiral phase transition or confinement. With this let- 
ter we go one step further including a transport equation 
self-consistently derived from the Nambu-Jona Lasinio 
(NJL) lagrangian. This allows to study microscopically 
the transport behavior of a fluid that includes the chiral 
phase transition looking at its impact on the relation be- 
tween the V2 and the rj/s of the system. The choice of the 
NJL is mainly driven by its wide and renowned applica- 
tion to study the QCD chiral phase transition by mean of 
effective lagrangians, even though the thermodynamical 
properties of QCD can be reproduced only qualitatively 
as briefly discussed in the following. 

The NJL Lagrangian is: 



C N jl = ipii^dn — mjip+g 



A 



a=l 



(1) 

with tp denoting a quark fields with Nf flavors ip = 
(u,d, ...)*, r a are the generators of the SU(Nf) group 
acting in flavor space with a = 1, (Nj — 1) . The 
to = diag(m u , rrid, ■•■) is the current Nf x Nf quark mass 
matrix in flavor space. In the following we will refer 
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to the Nf = 2, N c — 3 for calculations. As well know 
the theory is non-renormalizable, hence a cut-off A has 
to be introduced as a free parameter. The numerical re- 
sults shown in the following are derived using the Buballa 
parametrization: A = 588 MeV, gA 2 = 2.88, m = 5.6 
MeV l 21| that among the variety of parameterizations 
entails a behavior of e, P, c 2 s closer to the 1QCD results. 

A transport theory for the NJL model has been derived 
in the closed-time-path formalism combined with the ef- 
fective action method |22j. The main steps of deriva- 
tion are to perform a Wigner transformation of the Dirac 
equation of motion and of the related gap-equation asso- 
ciated to the Cnjl, Eq.([T]). Then one exploits the semi- 
classical approximation widely used for applications in 
heavy-ion collisions [23|, HH evaluating the expectation 
value of the four-point fermion interaction in the Hartree 
approximation (i.e. at mean field level). Finally only 
the scalar and vector components of the Wigner func- 
tion are retained thanks to the spin saturated nature of 
the systems we are interested in. One finally obtains 
the Boltzmann-Vlasov transport equations for the (anti- 
) quark phase-space distribution function / ± : 



pfd^fap) + M(x)d^M(x)d^f ± (x,p) = C(x,p) (2) 

where C(x,p) is the Boltzmann-like collision integral, 
main ingredient of the several cascade codes already de- 
veloped [III [2(|. We notice that respect to the al- 
ready implemented cascades the NJL dynamics intro- 
duces a new term associated to the mass generation. 
Also Eq.Q is formally the same as the widely used rel- 
ativistic transport approaches for hadronic matter like 
RBUU,uRQMD, RLV HI El Hj| , but with a vanishing 
vector field. However a key difference is that particles 
do not have a fixed mass and a self-consistent deriva- 
tion couples Eq. ^ to the mass gap equation of the NJL 
model that extended to the case of non-equilibrium can 
be written as: 



M(x) — m 
4<?A C 



= M(a 



d 3 p 1-f (x,p)- f+(x,p) 



(2tt) 3 



E p (x) 



(3) 



and determines the local mass M{x) at the space-time 
point x in terms of the distribution functions f ± (x,p). 

Eqs.([2]) and © form a closed system of equations 
constituting the Boltzmann-Vlasov equation associated 
to the NJL Lagrangian that allows to obtain self- 
consistently the local effective mass M{x) affecting the 
time evolution of the distribution function f ± (x,p). A 
seminal work on the transport equation associated to the 
NJL dynamics was done in Refs.[29|, [3(J, but without a 
collision term, not at finite r\j s and never applied to the 
physical conditions of ultra-relativistic heavy-ion colli- 
sions. 

For the numerical solutions of Eqs.© and ([3]) we 
use a three dimensional lattice that discretize the space 
as described in Ref.fllL HH- The standard test par- 
ticle methods that sample the distribution function / 



by mean of an ensemble of points in the phase-space 
is employed. The normalization condition is given by: 
JdTf ± =uj A(A) = N q (Nq) with T the phase space, A a 
(A a ) the number of test particles (antiparticles) which 
are inside the considered cell and uj the proper normal- 
ization factor that relates the test particles to the real 
particle number. 

In such a way it is possible to get a solution of the 
transport equations propagating the momenta of the test 
(anti-)particles by mean of the relativistic Hamilton's 
equation. For the numerical implementation they can 
be written in the discretized form as: 



M (r t) - 

P l {t + ) = Pi(t-) -2 5t "} 1 V r M a (n, t) + coll. 



r;(t+) = n(t-) + 25t 



Ei{t) 

Pi(t) 
Ei(t) 



(4) 



with t^ = t ± 5t and St the numerical mesh time. The 
term coll. on the right hand side of Eq.Q indicates the 
effects of the collision integral as described in Ref . (ill |26| . 
By mean of a reiterating procedure on time steps one gets 
the solutions of the transport equation coupling Eqs.(j4]) 
with the gap equation, Eq.© that discretized on lattice 
and for point-like test particles becomes: 



M„ 



8gN c 



M 



d 3 p 1 uj 
(2tt) 3 E ~ AT4 



(5) 



where A V a — tAt tanh r\ a is the volume of each cell 
of the space lattice given by At — 0.5/m 2 the area in 
the transverse direction and r] a the space-time rapidity 
of the center of the cell. The integral is instead the vac- 
uum contribution to the gap-equation which is a diver- 
gent quantity and it is regularized by a cutoff, A and has 
a simple analytical expression. 

The space-time dependence of the mass M a (r 7 t) = 
rn — 2(7^-0^^ influences the momenta of the particles be- 
cause the finite gradient of the condensate generates a 
force which changes the momentum of a particle pro- 
portionally to V r (V5V)j see Eq.Q. The last is negative 
because the phase transition occurs earlier in the sur- 
face of the expanding QGP fireball. Therefore the phase 
transition which take place locally results in a negative 
contribution to the particle momenta that makes the sys- 
tem more sticky respect to a free massless gas. 

Shear viscosity to entropy density ratio - The effect of 
the NJL mean-field can be evaluated looking at the so- 
called interaction measure normalized by T 4 , e ~4 P , that 
gives the deviation from the free gas relation between 
the energy density e and the pressure P that is also a 
measure of the breakdown of conformal symmetry. In 
lattice QCD it is known to be quite large with a peak 
at T c and a non negligible value up to T ~ 2 — 3T C . 
In Fig[T] (left) its behavior is shown for three different 
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FIG. 1: Left: Interaction measure shown as a function of 
temperature T for three different NJL parameter sets. Right: 
The rj/s and (,/s for the Buballa parametrization as a function 
of T for various cases as indicated in the legend. 



NJL parameterizations. Quantitatively the QCD behav- 
ior is correctly reproduced. In the following we will use 
the Buballa parametrization that gives the largest values 
closer to 1QCD where however the peak is found to be 
about a factor two larger at T c . This means that the 
numerical results on the impact of the mean field in the 
dynamics of RHIC collisions is reduced respect to a more 
realistic case. 

Our final goal is to study a fluid at fixed shear viscosity 
to entropy ratio 77/s extending the study started at cas- 
cade level [l(3, [H| • There the strategy was to normalize 
locally the cross section in order to fix the 77/s according 
to the simple relation a ■ 77/s =< p > /15n valid for a 
massless gas. Here because of the NJL field the particles 
acquire a mass hence both the viscosity and the entropy 
density are modified respect to the simple massless case. 
We briefly discuss the rj and s for a system of massive 
particles deriving the pertinent formula to renormalizc 
the cross section a in order to keep fixed locally the rj/s. 
Both 77 and s have been derived for a thermodynamical 
system and has been studied also for the case of the NJL 
model [3l[ . Here we derive expressions in terms of quanti- 
ties that can be used easily also in the numerical solution 
of the transport equations. A widely used formula for 77 
is deduced from the relaxation time approximation, like 
in Ref.|3l]. After integration by parts it is possible to 
write the shear viscosity for the general case of massive 
relativistic particles in terms of average quantities that 
can be easily evaluated numerically: 
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E 



\e 3 



(6) 



where r = \n(atrV r ei)\ is the relaxation time, i.e. 
the time interval between two collisions, and n is the 
total local density. One can easily see that in the ultra- 
relativistic limit (M — > 0) the well known formula for the 
shear viscosity, 77 = j%(p)/crtri is recovered. 

The entropy density cannot be related simply to the 
local density, s ~ An, as for the massless case. A suit- 
able method to evaluate locally s during the dynamical 



evolution of the collision is based on the use of the ther- 
modynamical relation sT = e + nT that shift the problem 
to the evaluation of e,n,T. These are easily calculated 
analytically (as for Fig[T]) but also numerically summing 
up the number of the test particles and their energy in 
each a— cell. To evaluate the temperature we exploit the 
general formula for a massive gas: 



e 

M 



Ki(z) 
K 2 {z) 



(7) 



with e the energy per particle, e/n, and K n (z) the 
modified Bessel functions of the second kind, z = 4f. 
We know e and M directly from the code hence we can 
use Eq.© to extract z = ^ and therefore the temper- 
ature. We have checked that the procedure works well 
performing calculation in a box with particles distributed 
according to a Boltzmann equilibrium distribution. 

The behavior of 77/s for a thermodynamical system is 
shown in Fig Q] (right) for a massless free gas (dashed line) 
and for the NJL (solid line). In both cases the cross sec- 
tion is fixed to a = 10 mb. The dot-dashed line indicates 
the lower bound for 77/s. We see that with a constant 
cross section a = 10 mb the 77/s is even lower than the 
1 /An at high T. The lower light lines show the behavior of 
bulk viscosity to entropy ratio £/s in the NJL model for 
two cases one with a = 10 mb (green solid line) as above 
and the other for an 77/s fixed at the 1/Air value (dot- 
dashed line). Of course for the massless case the bulk 
viscosity is zero, while for non- vanishing masses there is 
a link between the 77 and C through the relaxation time 
t. The last results in a smaller growth of £/s when 77/s 
is fixed respect to the case when cr tr is fixed. However in 
both cases we can see that only at T < 1 . 1 T c we have a 
non- vanishing £/s. This is due to the fact that the £ is 
expected to be proportional to the deviation of the sound 
velocity from 1/3, (c 2 s — 1/3) 2 that in the NJL model is 
known to occur only very close to T c . More importantly 
for our purposes is that in the NJL the C/s remains order 
of magnitudes lower that first extrapolation from 1QGD 
[32| and also much lower than the smaller values used for 
first studies with viscous hydrodynamics [19J. Consider- 
ing that even for much larger values of £/s hydrodynam- 
ics show a small effect, we can judge safe to discard any 
role of C/s in the following results for the elliptic flow. 

We notice that Eqs.© and ([7|) supply the formula for 
the normalization of the cross section in each ev— cell in 
order to keep fixed 77/s of the system: 



fir, a 



1 T a A(p 2 /E) a + M*(pyE 3 ) c 
15 (v re i) (e a + n a T a ) 77/s 



(8) 



for a massless gas M — > (p = E, e ~ 3nT) and Eq.© 
reduce to the simple relation atr,a = jf < P >a /n a 
used in Ref.sflFj, [ill, 0, HI for 77/s = 1/4tt. Eq.© will 
allow to extend such studies of a fluid at finite viscosity 
to the case of partons with finite mass. 

Elliptic Flow - We have run the simulations for Au+Au 
at y/~s~N~jy = 200 AGeV and 6 = 7 fm. The density dis- 
tribution in coordinate space is given by the standard 
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Glauber model. The maximum initial temperature is 
T = 340 MeV and the initial time is r = 0.6 fm/c as 
usually done also in hydrodynamical calculations. We fol- 
low the dynamical evolution of quarks, anti-quarks and 
gluons. The last has been included, even if they are not 
explicitly present in the NJL model, with the aim of using 
a realistic density for both the total and the (anti-)quark 
density in the simulation of the collisions. However glu- 
ons do not actively participate in the evaluation of the 
chiral phase transition, but they simply acquire the mass 
of the quarks not contributing to its determination ac- 
cording to the the NJL model. The justification for this 
choice relies on the quasi-particle models that are fitted 
to 1QCD thermodynamics [HI, One finds a similar 
behavior of M(T) for both gluons and quarks approxi- 
mately. Of course for a more quantitative calculation a 
more careful treatment would be needed but it is not rel- 
evant to the main objective of the present seminal work, 
considering also that anyway the NJL model cannot be 
used for an accurate quantitative study. 

In Fig. [5] (left) it is shown the time evolution of the 
average elliptic flow (v 2 ) for a constant transport cross 
section of cr tr =10 mb a typical value that is able to re- 
produce the amount of v 2 observed in experiments [2fj| . 
Comparing the two solid lines (black and green) we can 
see that the NJL mean field cause a decrease of (V2) of 
about 15%. The reduction of (V2) can be expected con- 
sidering that the NJL field produce a scalar attractive 
field that at the phase transition results in a gas of mas- 
sive particles. In Fig(3J we show also the elliptic flow 
at freeze-out as a function of the transverse momentum 
pr- One can see that the role of the mean field even 
increases with momentum affecting also particles at a 
Pt quite larger than the energy scale of the scalar con- 
densate (ipijj) ~ 300MeV. This is due to the fact that 
a high-pT particle collides mainly with the much more 
abundant particles in the bulk. These have an aver- 
age momentum comparable to the strength of the scalar 
field: pr ~ 2T ~ M c . Therefore the effect of the scalar 
field extends thanks to collisions into a range quite larger 
than one would naively think and the interplay between 
collisions and mean field is fundamental. We find that 
the presence of an NJL-field that drives the chiral phase 
transition suppress the V2(pt) by about 20% at pr > 1 
GeV. This would imply the need of a parton scattering 
cross section atr even larger than that estimated with 
the cascade model which was already quite larger than 
the pQCD estimates [5^, [36J. On the other hand the 
mean field modifies both the local entropy density re- 
duced by the mass generation, and the shear viscosity 
that increases as shown in FigQ] (right). Therefore even 
if the cross section is the same with and without the NJL 
field the system evolves with a different rj/s. Considering 
that one of the main goal is to determine the rj/s of the 
QGP we have investigated what is the action of the mean 
field once the rj/s of the system is fixed to be the same 
by mean of the cross section renormalization according 
to Eq. (|5]l. Therefore we have run simulations with and 



without the NJL-field but keeping constant locally the 

n/s. 
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FIG. 2: Left: Average elliptic flow as a function of time for 
Au + Au collisions in the mid- rapidity region \y\ < 1 at b= 7 
fm. Right: time evolution of the average transverse momen- 
tum for the same case as in the left panel. 
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FIG. 3: As in Fig. 2 for the elliptic flow as a function of pr- 

The results for Airrj/s — 1 are shown by dashed lines 
in Fig.s[2]and[3l the (black) dashed line is the case with 
only the collision term (cascade) while the (green) dashed 
line is the case with the field. We can see that once the 
rj/s is fixed there is essentially no difference in the calcu- 
lations with and without a field dynamics included. This 
is seen for both the average (vz) and the V2(pt) in Fig.s[5] 
(left) and [3] respectively. This is a key result that shows 
that even in a microscopic approach that distinguishes 
between the mean field and the collisional dynamics the 
V2(pt) is mainly driven by the rj/s of the fluid. In other 
words we have found that in a microscopic approach the 
rj/s is the pertinent parameter and the language of vis- 
cous hydrodynamics is appropriate. Of course this does 
not mean that V2{pt) in the transport theory is the same 
of the viscous hydrodynamics one, but that, and even 
more importantly, the direct relation between t>2 (pt) and 
rj/s is quite general and their relation is not modified by 
the NJL field dynamics. We have checked that this is 
valid also at other impact parameter (6 = 3, 5, 9 fm) and 
for larger rj/s up to tt~ 1 . This is of course very impor- 
tant for the determination of rj/s by mean of the data 
on elliptic flow and confirms the validity of the studies 
pursued till now even if they miss an explicit mean field 
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dynamics and/or the chiral phase transition. 

However we notice that while the i^Cpt) appears to 
be totally independent on the presence of the NJL-field 
once the rj/s is kept fixed, the time evolution of (1)2) 
still shows a slightly reduced elliptic flow at t ~ 3 — 
5/m/c. A similar difference can be observed also in the 
time evolution of the transverse momentum (px) shown 
in Fig. [2] (right). One may ask what is the physical origin 
of such differences. In principle there are two parameters 
affecting the i>2 (.Pt) and the (pt) : the sound velocity c s 
and the bulk viscosity £. As discussed previously it is 
safe to discard the possibility of any significant influence 
of the finite (/son our results considering its tiny value in 
the NJL, see FigJT] (right). It is instead reasonable that 
the weak decrease is due to the decrease of the sound 
velocity for NJL at T < 1.1 T c . It is well known that c 2 s 
decreases from 1/3 that is the value of a massless free 
gas and that this cause already in ideal hydrodynamics 
a decrease of the elliptic flow (37j |. On the other hand 
when the bulk of the system reaches this region most of 
the v 2 {pt) has already been built-up hence the effect of a 
moderate decrease of c s is quite weak and could explain 
the small difference still visible in the time evolution of 
(v 2 ) and (p T ). 

Conclusion - The novelty of the present work is to be 
the first study within a transport approach of a fluid at 
finite 77/s that includes the field dynamics of the chiral 



phase transition. Generally we find that at fixed cross 
section the effect of the NJL field is to reduce the elliptic 
flow by about a 20%. More importantly we can state 
that the presence of the NJL dynamics does not change 
the relation between the elliptic flow and the 77/s that 
remains the same as in the cascade models and at low pt 
is very close to the one from hydrodynamics (9I-H1I [33| . If 
such a finding is confirmed also for a more general class of 
interacting quasi-particle models it will make much safer 
and solid the determination of 77/s by i>2(pt)- In fact as 
we have shown the relation is independent on the micro- 
scopic details of the interaction once the EoS and/or the 
Cg(T) has been fixed. This will be investigated in the 
next future, in fact the kinetic theory and the numerical 
implementation presented here can be easily extended to 
quasi-particle models that are fitted to reproduce the en- 
ergy density and pressure of 1QCD results. In such a case 
it will be possible also to study the elliptic flow with a 
realistic behavior of c s (T) and the effect of a finite and 
sizable £/s on the elliptic flow complementing the study 
from viscous hydrodynamics that are subject to problems 
for not too small £/s and/or for pt > 3T 19]. 
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